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We show that a system of hard discs confined to a narrow channel exhibits a fragile-strong fluid 
crossover located at the maximum of the isobaric heat capacity and that the relaxation times for 
different channel widths fall onto a single master curve when rescaled by the relaxation times and 
temperatures of the crossover. Calculations of the configurational entropy and the inherent structure 
equation of state find that the crossover is related to properties of the jamming landscape for the 
model but that the Adam-Gibbs relation does not predict the relaxation behavior. We also show 
that a facilitated dynamics description of the system, where kinetically excited regions are identified 
with local packing arrangements of the discs, successfully describes the fragile-strong crossover. 



Upon cooling or compression, many materials, includ- 
ing supercooled liquids[l|, 0, gels and polymers, form 
amorphous, glassy solids, where the time needed for the 
system to structurally rearrange becomes longer than the 
experimental measurement. However, despite the ubiqui- 
tous appearance of glasses in nature and science, a com- 
prehensive understanding of this glass transition remains 
elusive. One useful approach has been to classify glass 
forming liquids on the basis of the temperature depen- 
dence of their structural relaxation timeSjT, on approach 
to the glass transition temperature T g j3j. In a strong 
liquid, t is linear in an Arrhenius plot of lnr vs 1/T, im- 
plying structural relaxation is a simple activated process. 
Fragile liquids exhibit super- Arrhenius behavior where r 
increases much faster and appears to diverge at positive 
temperatures. This suggests relaxation is cooperative 
and has led to speculation that there is a thermodynamic 
ideal glass transition underlying the kinetic behavior ob- 
served in experiment and simulation . 

Silica @, i| , silicon Q and water [ 



111 ] all appear to 



go through a dynamical crossover from a fragile liquid 
to a strong liquid at a crossover temperature, T x , that 
coincides with the Widom line, marked by a maximum 
in the heat capacity. However, a recent analysis of the 
transport coefficients of 84 different glass formers shows 
that the fragile-strong (FS) crossover occurs more widely 
than originally thought and suggests that the crossover 
temperature may be more relevant to the general features 
of dynamical arrest fl2| than T g , which is based on an 
arbitrarily chosen experimental relaxation time. 

In this letter, we study two contrasting paradigms used 
to describe glassy behavior: The first is the inherent 
structure landscape (ISL) [l3, 14 1, or the density packing 
landscape, which is the hard particle equivalent to the po- 
tential energy landscape [lj| , and the second is the facili- 
tated dynamics dUQjl (FD) approach. Both approaches 
are shown to capture key elements of the crossover but 
we also find that the local packing arrangements of discs, 
which ultimately give rise to the ISL, can be related to the 



kinetically excited regions appearing in the FD paradigm, 
providing a connection between the two approaches. 

Our model consists of N two-dimensional (2D) hard 
discs, with diameter a, confined between two hard walls 
(lines) of length L separated by a distance 1 < Hd/cr < 
1 + \/3/4. The particle-particle and particle-wall inter- 
action potentials are given by 
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respectively, where = |rj — i"i| is the distance between 
particles, r y is the component of the position vector for a 
particle perpendicular to the wall and ho = Hd — a. The 
occupied volume is cf> = Nira 2 / (ALHd). 

The exact partition function (l8| and complete jam- 
ming phase diagram (l9l . |20| for the system are known, 
making it an ideal tool for exploring the relationships 
between thermodynamics, dynamics and the ISL. Fig. Q] 
shows the four locally jammed packing configurations 
of the discs that can be combined to give the col- 
lectively jammed [2~lj ] inherent structures. The occu- 
pied volume fraction of the jammed states is <pj = 



7T/{ 



±H A 



(1 - 6) s/{2 - H d ) H d }, where is the 

mole fraction of defects (type 2 and 4 bonds). The 
most dense jammed state, 4>Jmax, occurs when = 0. 
The least dense jammed packing occurs when = 0.5 
and consists of a repeating unit cell of— 1 — 4 — 3 — 2— 
bonds because placing two defect bonds together (-2-2- 
, or -4-4-) results in an unstable configuration. Recent 
studies of this model found that configurations of the 
ideal gas mapped to jammed states, 4>ji g , correspond- 
ing to the maximum in the packing distribution with 
9 ig = 1/2 - 5^ 2 /10. The equilibrium fluid then only 
samples basins with a higher <f>j as (j> is increased fl9H22j]. 

We obtain the exact equation of state (EOS) using the 
transfer matrix method developed by Kofke et al [l8| . If 
the positions of the discs are fixed in the y-direction, the 
configurational integral in the cc-direction can be treated 
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FIG. 1. Local packing arrangements that can be combined 
to give collectively jammed states. Dashed lines connect the 
centers of neighboring discs in contact and the numbers iden- 
tify different "bonds". Bonds 1 and 3 are the locally most 
dense states. Bonds 2 and 4 represent the defect states. 



as a ID mixture of hard rods on a line. Taking the 
Laplace transform, gives the partition function in the 
N, P, T ensemble as 



A(N,P,T) 



A DN (f3P) N + 1 



dyK N {y,y). (2) 



Here, A is the thermal wavelength, P is the longitudinal 
pressure and K(yi,y 2 ) = exp[-PhoL x (y-L, y 2 )], with yi 
and 2/2 being the y-coordinates of two adjacent discs in 
contact. L x is the projection of the distance between the 
two contacting discs along the x-axis and is a function 
°f Hi j 2/2 and ho ■ Solving the eigenvalue problem asso- 
ciated with Eq. [2] yields the equation of state for the 
fluid in the thermodynamic limit. This gives us access 
to the heat capacity given by C p /Nk = (dH/dT) p = 
1 + Z/(1 + din {Z} /din {(j)}), where H = NkT + PV is 
the enthalpy in the hard discs system, Z = PV/NkT is 
the compressibility and V = h^L is the volume accessi- 
ble to the disc centers. Fig. [2ja) shows that C p /Nk as 
a function of T has a maximum. However, the C p maxi- 
mum in this system results from the binomial density of 
states and is not connected to the presence of an under- 
lying critical point associated with a liquid-liquid phase 
transition. 

Adam and Gibbs [23| argued that the rapid slow down 
in the dynamics of fragile liquids results from the decrease 
in accessible configurations at low temperatures or high 
densities and predicted relaxation times to behave as, r = 
Aexp (B/TS C ), where A and B are effectively constant, 
S c = kln(Nj) is the configurational entropy and Nj is 
the number of inherent structure basins accessible to the 
equilibrium fluid. The Adam-Gibbs relation predicts a 
divergence in r as Nj — > 1, causing S c — > 0, and it has 
been used to describe the relaxation in a wide variety of 
materials 

In our model, the configurational entropy fl9| is given 
by S c (<f>) /Nk = (1 - 9) ln(l - 9) - In 6 - (1 - 20) ln(l - 
26), where the equilibrium value of 6(<p) can be ob- 
tained by analytically quenching the fluid to its lo- 
cal inherent structure using the transfer matrix method 
and the information about the local geometry of four 
discs in contact contained in the chain product matrix 
K(y t ,y m )K(y m ,y n )K(y n ,y 3 ) Fig. G^a) plots S c as 
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FIG. 2. Cp/Nk for different values of Ha as a function of (a) 
(cPPV/NkT)- 1 and (b) 9. 




FIG. 3. (a) The equilibrium S c (b) Pis and (c) P v n as a 
function of 4> for different values of Hd- 



a function of <fi and shows that the rate of configurational 
entropy loss increases with increasing 4> at low densities 
but the impending Kauzmann catastrophe, thought to 
occur in fragile liquids, is avoided when S c goes through 
an inflection point. Consequently, the fluid has no ideal 
glass transition and S c only goes to zero in the limit 
4> <P.jmax and PV/NkT — > oo. This is the expected 
behavior for a system where the distribution of packings 
is determined by localized point defects [27} . In addition, 
Fig. 0(b) shows that the C p maxima, for the different Hd, 
all occur at the same 9 = 0.044 ± 0.002, which suggests 
it is the number of defects that controls the thermody- 
namics associated with the C p maximum. 

To calculate the relaxation times for the system as a 
function of <f>, we use molecular dynamics (MD) Simula- 
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tions where the time, t, has units of a^m/fcT) 1 / 2 . At any 
time during a MD trajectory, an instantaneous configu- 
ration can be mapped to its local inherent structure by 
considering the position of each disc relative to its two 
neighbors and using the triangular constraint, involving 
three neighboring discs, introduced by Speedy [24j. If the 
central disc lies below the line joining the centers of the 
outside discs, the central disc will pack at the bottom 
of the channel, otherwise it will pack at the top. This 
is equivalent to the analytical quench used to obtain S c . 
Once the local packing positions have been identified, the 
bonds between neighboring discs are assigned their la- 
bels, 1-4 (see Fig. CQ). The fluid remains within the basin 
of a single inherent structure for a short time before a lo- 
cal rearrangement of the discs, which changes the identity 
of some of the bonds, moves the system to a new inherent 
structure. We measure R(t), the fraction of bonds that 
have not changed at least once in time t as a function of 
t, and define the relaxation time as r = J Q R(t)dt. 

Our simulations were performed using N = 2000 discs 
and periodic boundaries in the longitudinal direction. 
At each <f> studied, 400V collisions were used to estab- 
lish equilibrium after the system had been compressed 
from the previous <j) using a modified version of the 
Lubachevsky and Stillinger algorithm (28| that maintains 
a fixed ratio of H^/a. Simulation lengths varied from 
200N collisions at low densities up to 10 6 A~ collisions at 
high densities and 80000 configurations were mapped to 
their inherent structure at each 0. R(t) always decays 
to zero in the time scale of the measurement, suggesting 
the system remains a fluid for all densities considered. 

For a hard particle system, <f)PV is a constant along 
an isobar and the Arrhenius law would predict that In r 
varies linearly with (j)PV/NkT. Fig. @] shows that r in- 
creases more rapidly than the Arrhenius law predicts at 
low densities (high T), which suggests the fluid is frag- 
ile, but we see a crossover to strong-fluid behavior at 
high densities, where the relaxation times increase lin- 
early. We also show fits of the data from the fragile region 
to the Vogel-Fulcher-Tammann [2<||3l[ (VFT) equation, 
r = Acxp [B/(T — Tvft)}, which predicts a divergence 
of the relaxation times at a temperature Tvft > OK, 
along with the parabolic law developed by Elmatad, 
Chandler and Garrahan 32|, l33j , which predicts no sin- 
gularity and is derived on the basis of the FD models. 
Both equations fit well when restricted to fragile fluid 
data and the Arrhenius equation provides the best fit for 
r above the crossover. Good fits of the VFT equation 
to a wide range of experimental and simulation data for 
supercooled liquids have been used as evidence for the 
presence of a thermodynamic singularity underlying the 
experimentally observed glass transition. However, we 
have already shown that our model does not exhibit an 
ideal glass transition, which suggests Tvft is simply a 
fit parameter with no physical significance. According to 
the Adam-Gibbs relation, lnr should vary linearly with 
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FIG. 4. Arrhenius plot for the relaxation times for H = 1.866. 
The dashed and dashed-dot lines represent fits to the data in 
the fragile region of the parabolic and VFT equations respec- 
tively. The solid line is the Arrhenius fit to the strong fluid 
region. Insert: The Adam-Gibbs plot for r. 



(f>PV/TSc, but the inset to Fig . |4] shows that this is not 
the case here. Sastry et al [34j recently found that the 
Adam-Gibbs relation did not hold in two-dimensions. 

The FS crossover occurs at the same </> as the maximum 
in the C p for all channel diameters, which is consistent 
with the studies that connect the crossover to the ther- 
modynamics of the Widom line. Furthermore, for each 
Hd we locate the temperature of the C p maximum, 

Tjnax j 

using our thermodynamic analysis, and define ro as the 
relaxation time at T max . The temperatures and relax- 
ation times are then rescaled by T max and ro respectively 
to give rise to the plot in Fig. 03a), where all the curves 
have collapsed onto a single master curve. Rescaling by 
any other temperature, for example, by defining an arbi- 
trary T g time scale, fails to collapse the data and leads to 
the impression that the systems with different Hd have 
different fragilities. This suggests T max = T x provides a 
more useful scaling temperature for our model. 

The relationship between the fragility of a fluid and 
the ISL has been explored in terms of the number and 
distribution of inherent structures 



351 but we can 



compare fragile and strong behaviors within the same 
model. In particular, we calculate the equivalent of the 
inherent structure pressure [III for a hard particle sys- 
tem directly from the configurational entropy as Pis = 
T(dS c /dL)u t h i then use the relation P = Pis + Pvib an d 
the exact EOS to isolate the vibrational contribution to 
the pressure (see Figs. [3jb) and[3]Jc)). Pjs makes a sig- 
nificant contribution to the overall pressure in the fragile 
fluid and reflects the fact that the configurational entropy 
of the system is varying rapidly as a function of density in 
the fragile fluid. However, Pjs then goes through a max- 
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FIG. 5. Master curve for the C p maximum rescaling of the 
relaxation times and temperatures for the Arrhenius plot (a) 
and the facilitated dynamics model (b). The solid lines rep- 
resent linear fits to the strong fluid region. 



imum at a density slightly higher than that of the C p 
maximum and rapidly decreases. The EOS of the strong 
fluid is entirely dominated by the vibrational component 
of the system rattling around the local jamming point of 
the inherent structures being sampled. This is consistent 
with experimental findings |37j that show strong behavior 
is dominated by density effects and local jamming while 
fragile relaxation is more thermally activated. 

The facilitated dynamics paradigm suggests relaxation 
and particle motion is driven by local microscopic dy- 
namical rules rather than any underlying thcrmodynam- 
A key ingredient of FD is the presence of ki- 



17] 



netically mobile regions that are able to influence the 
motion of neighboring regions, leading to the formation 
of chains of mobile particles in space-time. In addition, 
the theory argues that directed particle motion plays an 
important role. If a kinetically mobile region can acti- 
vate or deactivate a neighboring region independent of 
any previous motion, it is considered to be directionally 
independent. Then the system behaves like a strong fluid 
and lnr w — hi Cfd, where Cfd is the concentration of 
kinetically excited regions. The expectation that struc- 
tural relaxation in a fragile fluid is cooperative is cap- 
tured by having directional correlation between the suc- 
cessive movement of particles in the kinetically excited 
regions. A FS crossover is predicted to occur when el- 
ements of both mechanisms are present in the system. 
While FD models have been parameterized to fit exper- 
imental data, only recently have there been efforts to 
identify the kinetically excited regions at a microscopic 



level [38H40J and most studies of FD have focused on spin 
models where the dynamic rules are included by construc- 
tion. 



In the current model, we are able to identify these ki- 
netically excited regions as the defects states (2, 4 bonds) 
in the jammed configurations. Relaxation in the system 
occurs in three ways: A) A particle next to an isolated 
defect can hop into the defect, causing the defect to move. 
This occurs with equal probability in both directions and 
leads to strong fluid behavior. B) Two neighboring de- 
fects moving toward each other create a local configura- 
tion with bonds —1 — 4 — 4 — 3— or —3 — 2 — 2 — 1— that 
is unjammed, leading to a spontaneous collapse of the 
central disc and an annihilation of the defects to form a 

— 1 — 3—1 — 3— locally jammed state. It is the spontaneity 
of the particle rearrangements associated with the defect 
annihilation, following the initial particle hop that brings 
the defects together, that is characteristic of the cooper- 
ative relaxation in a fragile fluid. C) A non-defect state 
can create two neighboring defects that move apart. The 
equilibrium number of kinetically excited regions in the 
system is then just Cfd = which we obtained from 
our analytical quench of the system. At low densities, 
there is a high concentration of defects that can interact 
and the directed creation-annihilation mechanisms dom- 
inate, giving rise to a fragile fluid that crosses over to 
a strong fluid as Cfd decreases below its critical value. 
Fig-EIb) shows that lnr/Vo becomes a linear function of 

— \xl(Cp° d /Cfd), where Cp° D is the number of defects in 
inherent structures of the ideal gas, 8 ig , highlighting the 
strong fluid behavior at low Cfd (high densities). The 
figure also shows that the relaxation times for this model 
collapse to a single master curve when rescaled by To, 
which, on the basis of Fig. EJb), is equivalent to rescal- 
ing by the relaxation time of the system containing the 
critical number of defects associated with the crossover. 

Our work suggests that the local packing environ- 
ments of particles, along with the way they interact, may 
serve as the important microscopic ingredients in the FD 
paradigm and points to a new analysis that can be ex- 
plored in higher dimensions. In the case of hard parti- 
cles, it may be useful to identify local packing structures, 
or local tilings, in the jammed inherent structures (4l| 
as defects. In more complex liquids, such as silica and 
water, the local geometry of a particle is dominated by 
the formation of a random tetrahedral network (RTN) 
of bonds. A recent study of the dynamics of ST2 wa- 
ter [l(| , both above and below T x , showed that the tem- 
perature dependence of the diffusion coefficient could be 
explained in terms of the concentration of local defects 
in the RTN. Similar results have been found in models of 
colloids [42} and nanoparticle systems 43 1 with highly di- 
rectional, tetrahedral bonding. It is still not known if the 
movement of particles in relation to the RTN defects can 
be described by the FD model, but these studies, along 
with ours, strongly suggest local packing and particle ge- 
ometries may play an important role in the dynamics of 
fluids in general. 
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